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Many features of granular media can be modeled by a fluid of hard spheres with inelastic colli- 
sions. Under rapid flow conditions, the macroscopic behavior of grains can be described through 
hydrodynamic equations accounting for dissipation among the interacting particles. A basis for 
the derivation of hydrodynamic equations and explicit expressions appearing in them is provided 
' by the Boltzmann kinetic theory conveniently modified to account for inelastic binary collisions. 

The goal of this review is to derive the hydrodynamic equations for a binary mixture of smooth 
inelastic hard spheres. A normal solution to the Boltzmann equation is obtained via the Chapman- 
Enskog method for states near the local homogeneous cooling state. The mass, heat, and momentum 
fluxes are obtained to first order in the spatial gradients of the hydrodynamic fields, and the set of 
^ . transport coefficients are determined in terms of the restitution coefficients and the ratios of mass, 

concentration, and particle sizes. As an example of their application, the dispersion relations for 
the hydrodynamic equations linearized about a homogeneous state are obtained and the conditions 
for stability are identified as functions of the wave vector, the dissipation, and the parameters of the 
^ ' mixture. The analysis shows that the homogeneous reference state is unstable to long enough wave- 

C/3 ' length perturbations and consequently becomes inhomogeneous for long times. The relationship of 

this instability to the validity of hydrodynamics is discussed. 

a 

^ ■ I. INTRODUCTION 

o 

^ ' Granular media under rapid flow conditions exhibit a great similarity to molecular fluids. This fact has stimulated 
the use of hydrodynamic- like type equations to describe the macroscopic behavior of such systems. The main difference 
from ordinary fluids is the absence of energy conservation, leading to both obvious and subtle modifications of the 
Navier-Stokes hydrodynamic equations. To isolate the effects of such collisional dissipation from other important 
1 properties of granular media, an idealized microscopic model system is usually considered: a system composed by 
T-j. ' smooth hard spheres with inelastic collisions. As in the elastic case, the collisions are specified in terms of the change 
[~««, . in relative velocity at contact but with a decrease in the magnitude of the normal component measured by a positive 
' restitution coefficient a < 1. This parameter distinguishes the ideal granular fluid (a < 1) from the ideal normal fluid 

■ Although many efforts have been devoted in the past few years in the understanding of granular fluids, the derivation 
of the form of the transport coefficients is still a subject of interest and controversy. In the low-density regime, 
this problem can be addressed by using the Boltzmann kinetic theory conveniently modified to account for inelastic 
collisions "2]. The conditions for hydrodynamics are expected to be similar to those for normal fiuids. For a given initial 
: state there are two stages of evolution. First, during the kinetic stage there is rapid velocity relaxation to a "universal" 
r-| ' velocity distribution that depends on the average local density, temperature, and flow velocity. Subsequently, the 
Q , hydrodynamic stage is described through a slower evolution of these local hydrodynamic fields as they approach 
O ' uniformity. The solution to the Boltzmann equation in this second stage is said to be "normal", where all space 
I and time dependence of the distribution function occurs through the hydrodynamic fields. The Chapman- Enskog 
• • method |3l| provides a means to construct explicitly the form of this normal solution as a perturbation expansion 
rN in the spatial gradients of the fields. This solution is then used to evaluate the fluxes in the macroscopic balance 
^ , equations in terms of these gradients. To lowest order the balance equations become the granular Euler equations; 
■ ■ ■ ' to second order they are the granular Navier-Stokes equations. In carrying out this analysis, explicit forms for the 
transport coefficients are obtained as functions of the restitution coefficient and other parameters of the collision 
operator. This derivation of hydrodynamics from the Boltzmann equation has been widely covered in the case of a 
monocomponent gas where the particles are of the same mass and size jlj. However, a real granular system is generally 
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characterized by some degrees of polidispersity in density and size, which leads to phenomena very often observed 
in nature and experiments, such as separation or segregation. Several attempts to apply the Boltzmann equation 
to derive transport coefficients for a multicomponent system began time ago |1| , but the technical difficulties of the 
analysis entailed approximations that limited their accuracy. In addition, all these works are based on the assumption 
of energy equipartition so that the partial temperatures Ti are made equal to the global granular teniperature T . 
Nevertheless, kinetic theory studies 0,0, 0,13, computer simulations P ITfl ITll| and even real experiments 0,0] have 
clearly shown the breakdown of energy equipartition. As a consequence, many of the previous results obtained for 
granular mixtures must be reexamined by using a kinetic theory which takes into account the effect of temperature 
differences on the transport coefficients. 

The goal of this paper is to derive the hydrodynamic equations of a granular binary mixture from the Boltzmann 
kinetic theory. These equations are derived by applying the Chapman-Enskog expansion about a local homogeneous 
cooling state (LHCS) that is analogous to the local equilibrium state for a gas with elastic collisions. In the first order 
of the expansion, the irreversible parts of the mass, heat, and momentum fluxes are calculated and the eight transport 
coefficients identified. These coefficients are expressed in terms of the solutions to a set of coupled linear integral 
equations fl^. The analysis carried out here is more complete than previous studies it is exact to leading order in 
the dissipation but not limited to weak dissipation. However, for practical purposes, the above integral equations are 
solved approximately by using the leading terms in a Sonine polynomial expansion. Such approach compares quite well 
with the results obtained from numerical solutions [^Il5|| of the Boltzmann equation by means of the Direct Simulation 
Monte Carlo (DSMC) method 16J. The explicit knowledge of the transport coefficients allows quantitative application 
of the nonlinear hydrodynamic equations to a number of interesting problems for mixtures, such as segregation and 
separation. Here, the simplest example of an application is considered, small perturbations of a spatially homogeneous 
state. The dispersion relations for the hydrodynamic modes are obtained from the hydrodynamic equations linearized 
about the homogeneous form of the LHCS. Linear stability analysis shows two shear modes and a heat mode to be 
unstable for long wavelength excitations. This means that small perturbations or fluctuations about the homogeneous 
state will grow. The implications of this for the derivation of the hydrodynamic equations are discussed here as well. 

The plan of the paper is as follows. In Sec. m we review the Boltzmann equation and associated macroscopic 
conservation laws. We also give a brief survey of the application of the Chapman-Enskog method around the LHCS. 
This state is analyzed in Section IlIII while in Sec. llVl we give the form of the Navier-Stokes hydrodynamic equations. 
Theoretical results are compared with simulation data for the temperature ratio and the shear viscosity coefficient. 
Section deals with an analysis of the stability of the linearized hydrodynamic equations. We close the paper in 
Section IVll with a discussion of the results presented. 



II. BOLTZMANN EQUATION AND CONSERVATION LAWS 

We consider a binary mixture of smooth hard spheres of masses mi and m2, and diameters ui and (72 • The 
inelasticity of collisions among all pairs is characterized by three independent constant coefficients of normal restitution 
dill oti2i and a\2 = a2i, where aij is the restitution coefficient for collisions between particles of species i and j. In 
the low-density regime, the distribution functions /i(r, v; t) (i = 1, 2) for the two species are determined from the set 
of nonlinear Boltzmann equations 

idt + VI . V) /. = 5] [vi|/,(t), j,[t)\ , (1) 
i 

where the Boltzmann collision operator Jy [vij/^, /j] describing the scattering of pairs of particles is 

Ji3\^-^\h.!'3\ = Cr^y" d-W-i J dCTe(S • gi2)(CT • gia) 

X [ar2/,(r,v'i,t)/,(r,vi,i)-/z(r,vi,t)/,(r,V2,i)] , (2) 

where aij = (cTj + (jj) /2, a is a unit vector along their line of centers, is the Heaviside step function, and gi2 = 
Vi — V2. The primes on the velocities denote the initial values {v'i,V2} that lead to {vi,V2} following a binary 
collision: 

v[ = vi - fiji {l + a^r.^) ■ gu)^ , V2 = V2-FAiij (l + a^j-^) (S-gi2)?, (3) 



where /i.y ~ rrii/ (rrii + rrij ) . 
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The relevant hydrodynamic fields are the number densities n^, the flow velocity u, and the "granular" temperature 
T. They are defined in terms of moments of the distributions fi as 

^ J dvf,{-v) , pu = Y^ J dvm,vf,{-v) , (4) 

i 

^T = P = J2j d^r^V'Mv), (5) 

where ti = rti +712 is the total number density, p = toi?ii +TO2't-2 is the total mass density, p is the hydrostatic pressure, 
and V = V — u is the peculiar velocity. At a kinetic level, it is convenient to introduce the kinetic temperatures Ti 
for each species defined as 

^n,r, = I d^^V^U (6) 

The collision operators conserve the particle number of each species and the total momentum but the total energy 
is not conserved: 

J d^J,, [v|/„ /,] - 0, Y.j f^vm,v J,, [v|/„ /,] = 0, (7) 

^ j dvim,yV.,[v|/„ /,] = -^nTC, (8) 

where C, is identified as the cooling rate due to inelastic collisions among all species. The macroscopic balance equations 
follow from the Boltzmann equation |^ and Eqs. Q and JSJ. They are given by 

Dtn, + niV • u + = , (9) 
nil 

Au + p"^VP = 0, (10) 

AT--^^ + A(v.q+P: Vu) = -CT. (11) 

i 

In the above equations, Dt — dt + u ■ \7 is the material derivative, 

j. = m,ydvV/,(v) (12) 
is the mass fiux for species i relative to the local flow, 

P = Y. f ^^vm,VV/,(v) (13) 

is the total pressure tensor, and 

/ rfvim,F2V/,(v) (14) 

i 

is the total heat flux. 

The utility of the balance equations (P|l- (|ll|l is limited without further specification of the fiuxes and the cooling 
rate, which in general have a complex dependence on space and time. However, for sufficiently large space and time 
scales, one expects that the system reaches a hydrodynamic regime in which all the space and time dependence is 
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given entirely through a functional dependence on the hydrodynamic fields n^, u, and T. The corresponding functional 
dependence of ji, P, and q on these fields are called constitutive equations and define the transport cocfhcients of 
the mixture. The primary feature of a hydrodynamic description is the reduction of the description from many 
microscopic degrees of freedom to a set of equations involving only five local fields. At a kinetic level, the constitutive 
equations are obtained when one admits the existence of a normal solution to the Boltzmann equation where the 
velocity distribution function depends on r and t only through its functional dependence on the fields, namely, 

/,(r,vi,t) = /4vi|ni(r,t),r(r,t),u(r,<)]. (15) 

The Chapman-Enskog method Q generates this normal solution explicitly by means of an expansion in gradients of 
the fields: 

+ ••• , (16) 

where e is a formal parameter measuring the nonuniformity of the system. The local reference state ff^ is chosen 
to give the same first moments as the exact distribution fi. The time derivatives of the fields are also expanded as 
dt — df^ + €d[^^ + • • • . The coefficients of the time derivative expansion are identified from the balance equations 
after expanding the fluxes, and the cooling rate ^ in a similar series as (|16|l . We close this section with 
some comments about the partial temperatures Ti introduced in Eq. JSJ). In contrast to the global temperature T 
these partial temperatures are not hydrodynamic fields and do not appear directly in the balance equations if^- lfTT)) 
and in the normal solution H15|) . However, as indicated in the next section, these partial temperatures characterize 
differences in the velocity distributions for the different species and therefore are important for their quantitative 
description. Ultimately, the transport coefficients obtained depend on the ratio T^jTi which is expressed in terms of 
the mechanical properties of the particle collisions, through the requirement that the solution be normal. 



III. LOCAL HOMOGENEOUS COOLING STATE 



To zeroth order in e, the kinetic equations become 
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The normal form of /j''^'' requires that all its time dependence occurs only n^, u, and T . Consequently, jf'^ must be 
of the form 

/f)(v)=n,t;o-3$,(U/«o) , (18) 

where v^it) = 2T(t)(m\ +r7i2)/ (rn\ra2) is a thermal velocity defined in terms of the temperature Tit) of the mixture. 
The balance equations to this order become d^f^ Xi = 9j°''u — OjT^^df^^T = — C^'^^. Here, Xi = rii/n is the 

mole fraction of species i and the cooling rate C^"^ is determined by Eq. ||SJ| to zeroth order. Hence, all the time 
dependence to this order is through the temperature and so, in dimensionless form, Eq. (|17|l becomes the following 
time independent equation for 

irg^-(vt<f.) = E4['J'-'f.]' (19) 

^ i 

where V* — V/vq, C* — C,i/'nai2VQ, and J,* — v^Jij /nnia'^2- Since the distribution functions are isotropic, it follows 
from Eqs. (|12|I - H13() that the zeroth order mass and heat fiuxes vanish while, for the same reason, the momentum flux 
is diagonal with a coefficient that is just the sum of the partial pressures, i.e., 

jr)=0, q(")=0, Pi?=p<5„,. (20) 

In summary, Eq. (|18|l provides the normal solution to lowest order in the gradients and use of Eq. H2U|I in the balance 
equations gives the corresponding Euler level hydrodynamics. Note that is determined self-consistently in the 
solution to Eq. . 

Just as for the one component gas case^J the exact form of $i has not yet been found, although a good 
approximation for thermal velocities can be obtained from an expansion in Sonine polynomials|5|. In the leading 
order, $i is given by 



3 /2 



(21) 



5 




FIG. 1: Plot of the coefBcients Ci versus a for a\/a2 = 1, = | and m\/m2 — 2. The solid line and the circles refer to ci 
while the dashed line and the squares correspond to C2. The dotted line and the triangles refer to the common value in the 
monocomponent case. The lines are the theoretical results and the symbols are the simulation results. 



Here, 9i = (jnij^i) ^ and 7^ — Ti/T. The coefficients Ci (which measure the deviation of $j from the reference 

MaxweUian) are determined consistently from the Boltzmann equation. An important observation is that the partial 
temperatures have been introduced here since they provide the actual mean kinetic energies of the distributions 
through Eq. ©. This would appear to contradict the assumption of a normal solution by introducing additional 
variables in the distribution functions. To see that this is not the case, the time evolution for the temperature ratio 
7 = T1/T2 is calculated to be 

5f ln7 = Cf -Cf^ (22) 
where we have introduced the cooling rates for the partial temperatures Ti as 

- / rfvim.yV,[v|/f ]. (23) 

The total cooling rate can be expressed in terms of the coolings cl"' as 

C(°)=r-iE^'^'Cf^- (24) 

i 

The fact that fi depends on time only through T(t) necessarily implies that the temperature ratio 7 must be 

independent of time and so, Eq. gives the condition (^f ^ — (^^^ = . In the elastic case, where °^ is a 
local Maxwellian, the above condition yields Ti{t) = T2{t) = T{t) and the energy equipartition applies. However, in 
the inelastic case, the equality of the cooling rates leads to different values for the partial temperatures, even if one 
considers the Maxwellian approximation to ff^'^ . Nevertheless, the constancy of 7 assures that the time dependence 
of the distribution is entirely through T(i), and in fact the condition of equal cooling rates assures that the partial 
temperatures can be expressed in terms of T. 

The explicit form of the approximation H21I) provides detailed predictions for the temperature ratio Ti /T2 (through 
calculation of the cooling rates) and for the cumulants Ci [from Eq. H19() ] as functions of the mass ratio, size ratio, 
composition and restitution coefficients0|. This completely fixes the approximate distribution function. In Fig. ^ 
we show the dependence of ci and C2 on the (common) restitution coefficient ay = a for uxja^ — 1, ^\ = ^, and 
mi/ni2 = 2. We also present the corresponding results obtained from a numerical solution |9| of the Boltzmann 
equation by means of the DSMC method|lS]. The agreement between the simulation data and the theoretical results 
is excellent. Further, the small values of the coefficients Ci supports the assumption of a low-order truncation in 
the polynomial expansion of the distribution functions. One of the main findings of our theory is that, except for 
mechanically equivalent particles, the partial temperatures are in general different (7 = T1/T2 ^ 1). This conclusion 
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FIG. 2: Plot of the temperature ratio 7 = T1/T2 versus a for 0^/02 = 1, a;i = | and three different values of the mass ratio: 
m\/m2 =0.1 (dotted line and triangles), mi/m2 = 2 (dashed line and squares), and m\/m2 = 10 (solid line and circles). The 
lines are the theoretical results and the symbols are the simulation results. 

contrasts with all previous resultsQ derived for granular mixtures, where it was implicitly assumed the equipartition of 
granular energy between both species (i.e., 7=1). The breakdown of the energy equipartition has been confirmed by 
computer simulations [gII^. 1 1 d 1 1 ll | and by real experiments of vibrated mixtures in two |l2j| and three |l3j | dimensions. 
For the sake of illustration, in Fig. |21we plot the temperature ratio 7 versus a for (Ti/(T2 = 1, xi = |, and several 
values of the mass ratio as given by the first Sonine approximation and by Monte Carlo simulations. The comparison 
between theory and simulation data shows an excellent agreement, even for strong dissipation. We also observe that 
for large differences in the mas ratio, the temperature differences are quite important. 

IV. NAVIER-STOKES HYDRODYNAMIC EQUATIONS 

Implementation of the Chapman-Enskog method to the first order in e has been carried out and without approx- 
imation recently Technical details of these calculations will not be repeated here and only the main results are 
offered. The solution of the Boltzmann equation to first order in the spatial gradients can be written as 



(25) 



The coefficients A, Bi, Ci, and "Di^M are functions of the peculiar velocity and the hydrodynamic fields. These 
coefficients obey inhomogeneous integral equations involving the linearized Boltzmann operators. Solubility conditions 
for the existence of solutions has been also proven. The corresponding constitutive equations found to this order are 



.(1) 
Jl 



mim2n 



DVxi - -Dp\/p 



Jl 



(26) 



,(1) 



-T^D"\/xi - L\/p - AVT, 



(27) 



p(i) _ _„ 



ViUk 



u 



(28) 



Furthermore, there is not contribution to the cooling rate at this order, i.e., C'-^-' = 0. This is a consequence of the 



symmetry of the distribution function /-[^■' since the latter does not contain any contribution proportional to the 
divergence of the flow velocity field. This property is special to the low density Boltzmann kinetic theory and such 
terms occur at higher densities The transport coefficients appearing in Eqs. H2t)|l - H28|) are the diffusion coefficient 
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D, the thermal diffusion coefficient D', the pressure diffusion coefficient Dp, the Dufour coefficient D", the thermal 
conductivity A, the pressure energy coefficient L, and the shear viscosity 77. These transport coefficients can be easily 
written in terms of 3i, Ci, and Vi^ke- As in the elastic caseQ, an accurate estimate of the transport coefficients can 
be obtained by using the leading terms in a Sonine polynomial expansion. Thus, the transport coefficients are known 
as explicit functions of the restitution coefficients and the parameters of the mixture (masses, sizes, and composition) 
and their expressions are not limited to weak inelasticity. 

These expressions for the mass flux, the pressure tensor, the heat flux, and the cooling rate provide the necessary 
constitutive equations to convert the balance equations into a closed set of six independent equations for 

the hydrodynamic fields. Since Eqs. (|26|) - (|28|l have been represented in terms of the gradients of the mole fraction 
Xi, the pressure p, the temperature T, and the flow velocity u, it is convenient to use these as the independent 
hydrodynamic variables. This means that, apart from the balance equations © and © for u and T, we need to get 
the corresponding balance equations for xi and p. These equations can be easily obtained from and and are 
given by 

Dtxi + ^ V-ji = 0, (29) 

Ap + pV-u+^(V-q+P: Vu) =-Cp- (30) 

Therefore, when the expressions of the fluxes H2t)|l - H28() . and the cooling rate C C^'^'' are substituted into the balance 
equations (jHl, Q, and lp!n| one gets a closed set of hydrodynamic equations for xi, u, T, and p. These are the 
Navier-Stokes equations for a mixture, which are given by 

D,x, = . ( "^^DVx, + P-D.Vp + ^D'Vt) , (31) 

{Dt + 0p+ f PV • u = ^ V • {T^D"\/xi + L\/p + AVT) 

+ ^?7 (^eUk + VfeU£ - ^4£V • V<>Ufc, (32) 

iD, + OT+'-pV-u^ :^DVx, + ^D,Vp + ^D'Vt) 

6 n \mi m2 / VP P ^ J 

+ J-V • (T^D"\/xi + L\/p + AVr) 

+ 77-^ ( V^Mfc + Vfc-u^ - |(5fe£V • u ) V^Mfe, (33) 
3n V 3 / 



Dtug + V^p = p ^S/kV yV£Uk + ^kUi - -5ki>^ ■ uj . (34) 

In general, the influence of dissipation on the transport coefficients is quite important even for moderate values of the 
restitution coefficient (say a ~ 0.9). Further, the effect of the additional contributions to the transport coefficients 
coming from the partial temperature differences can be quite significant, even for weak dissipation. These new 
contributions were not taken into account in previous work on granular mixtu res M . The accuracy of the approximate 
Sonine solution to the resulting integral equations has been recently tested [1^ at the level of the shear viscosity 
coefficient. Specifically, the DSMC method has been used to solve the Boltzmann equation in the special hydrodynamic 
state of uniform shear flow. An external thermostat is introduced to compensate for the energy lost in collisions so 
that the shearing work still heats the system. In the long time limit, the Navier-Stokes shear viscosity coefficient can 
be measured from simulations and compared with the predictions given by the Chapman-Enskog method. In Fig. |31 
we plot the ratio r]* {a) / rj* {1) versus the mass ratio mi/m2 for <J\Ig2 = 1, Xi = i, and three different values of the 
(common) restitution coefficient a. Here, 77* = mj/riT, with 1/ — y^ncr^g^o and r?*(l) refers to the elastic value of the 
shear viscosity coefficient. Again, the symbols represent the simulation data while the lines refer to the theoretical 
results obtained from the Boltzmann equation in the first Sonine approximation. We see that in general the deviation 
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FIG. 3: Plot of the ratio r]* (a) / rj* (1) as a function of the mass ratio mi/m2 for axjai — ni/n2 ~ 1 and three different values 
of the restitution coefficient a: a = 0.9 (circles), a = 0.8 (squares), and a = 0.7 (triangles). The lines are the theoretical 
predictions and the symbols refer to the results obtained from the DSMC method. 




FIG. 4: Plot of the ratio tj* {a)/ri* (1) as a function of the size ratio cri/(T2 for m\/m2 = 4, n\/n2 = 1 and three different values 
of the restitution coefficient a: q = 0.9 (circles), a = 0.8 (squares), and a — 0.7 (triangles). The lines are the theoretical 
predictions and the symbols refer to the results obtained from the DSMC method. 

of ri*{a) from its functional form for elastic collisions is quite important. This tendency becomes more significant 
as the mass disparity increases. The agreement between the first Sonine approximation and simulation is seen to be 
in general excellent. This agreement is similar to the one previously found in the monocomponent case '2^,'2T1. At 
a quantitative level, the discrepancies between theory and simulation tend to increase as the restitution coefficient 
decreases, although these differences are quite small (say, for instance, around 2% at a = 0.7 in the disparate mass 
case TOi/m2 = 10). The influence of the size ratio on the shear viscosity is shown in Fig. ^ for mi/m2 = 4 and 
xi = i. We observe again a strong dependence of the shear viscosity on dissipation. However, for a given value of 
a, the influence of ai/a2 on rj* is weaker than the one found before in Fig. |21for the mass ratio. The agreement for 
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both a = 0.9 and a = 0.8 is quite good, except for the largest size ratio at a = 0.8. These discrepancies become 
more significant as the dissipation increases (say a = 0.7), especially for mixtures of particles of very different sizes. 
In summary, according to the comparison carried out in Figs. El and 01 one can conclude that the agreement between 
theory and simulation extends over a wide range values of the restitution coefficient, indicating the reliability of the 
first Sonine approximation for describing granular flows beyond the quasielastic limit. 



V. LINEARIZED HYDRODYNAMIC EQUATIONS AND STABILITY 



The simplest application of the hydrodynamic equations is to small initial perturbations about the strictly homo- 
geneous state for a large system. Then the nonlinear equations can be expanded to linear order in the deviations of 
Xi, u, T, and p about their homogeneous values. This leads to partial differential equations with coefficients that 
are independent of space but they depend on time since the reference state is cooling. As in the monocomponent 
caseQ, this time dependence can be eliminated through a change in the time and space variables, and a scaling of 
the hydrodynamic fields. Let Sya{r,t) = ya(r,t) — yHa{i) denote the deviation of {xi,u, T,p} from their values in 
the homogeneous state. We introduce the following dimensionless space and time variables: 







(35) 



where vonit) and vnit) has been introduced before. The superscript H indicates that the quantity is evaluated in 
the homogeneous cooling state. This means that 



dtXiH = 0, uh = 0, dt InTff = dt \npH = -(h 
A set of Fourier transformed dimensionless variables are then defined by 

3-*-^,52/„(^,r), 



(36) 



(37) 



, . Sxikir) (5uk r) 

XlH VoH[T) 



(38) 



Ph{t) ' 



(39) 



In terms of the above variables the linearized hydrodynamic equations for the set {pk, Wk, ^k, Hk} separate into a 
degenerate pair of equations for the transverse velocity components Wki (orthogonal to k) 



d C 



ri*k^ I Wk 



0, 



2 

and a coupled set of equations for pkj ^k: Ilki and the longitudinal velocity component z«k|| (parallel to k) 



dr 



(40) 



(41) 



where now Syj^aiT) denotes the four variables (pk, ^k, Hk, Wkll) • The matrices in this equation are 













-C* 



(42) 



^/ 



/O 

















i^O 


— M21 
Xifl-\-X 



(43) 



10 



-D* ~x-^d'* -x^^d; 

-xi(^D"* h^D*] ^-f D* ~ |A* + Dl 

y3 xifi-\-X2 J xifi-\-X2 J 3 3 xifi-\-X2 P 

-^xiD"* -|A* -|L* 

V -|77*; 



(44) 



3' 

In these equations, fi = mi/m2, Xi = riiH /nu, Q* ~ C'h' I^h- Moreover, we have introduced the reduced Navier-Stokes 
transport coefhcients 



PhvIh ' 



(46) 



i^"-=-^^%^: L*.^. r^^. (47) 

Equation H40|l is decoupled from the rest and can be integrated directly leading to 

Wk±(T") = Wk±(0)exp(sir), (48) 

where 

s± = IC - Vk\ (49) 

This identifies two shear modes. The remaining modes have the form exp(s„T) for n = 1,2,3,4, where s„ are the 
solutions of the corresponding quartic equation. A detailed study on the dependence of the hydrodynamic modes 
on the inelasticity and the parameters of the mixture will be given elsewhere. The linear hydrodynamic equations 
characterize the stability of the homogeneous cooling state. We see from Eq. (I49|l that s± for the two shear modes is 
positive for A; < fc^, where 



^* \ 1/2 



(50) 



Thus, initial long wavelength perturbations of the HCS that for instance, excite the shear mode will grow exponentially, 
representing an instability of the reference state. 

The wave vector dependence of the remaining four modes is more complex. However, the long wavelength stability 
can be obtained from the eigenvalues of M'^°) which are easily found to be 

4c*,o,o,ir)- (51) 

Hence, there is another unstable mode. In summary, at asymptotically long wavelengths the spectrum of the linearized 
hydrodynamic equations is comprised of a decaying mode at — C*/2, a two-fold degenerate mode at 0, and a three-fold 
degenerate unstable mode at C*/2. 

VI. DISCUSSION 

The primary objective of this review has been to obtain the hydrodynamic description of a binary mixture of 
granular gases from the Boltzmann kinetic theory. The Chapman-Enskog method is used to solve the Boltzmann 
equation up to the Navier-Stokes order and the associated transport coefficients are given in terms of the solutions of 
a set of linear integral equations. A practical evaluation of these coefficients is possible by taking the leading terms in 
a Sonine polynomial expansion. In contrast to previous studies Q , our results apply for arbitrary degree of dissipation 
and they are not restricted to specific values of the parameters of the mixture. The explicit knowledge of the transport 
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coefHcients and the cooling rate allows one to perform an study of the linearized hydrodynamic equations around the 
homogenous cooling state. The stability analysis shows that the homogeneous cooling state is unstable to long enough 
wavelength perturbations. 

The reference state in our Chapman-Enskog expansion has been taken to be an exact solution of the uniform 
Boltzmann equation. An interesting and important result of this solution is that the partial temperatures of each 
species must be different. This does not mean that there are additional hydrodynamic degrees of freedom since their 
cooling rates are equal and consequently, the partial temperatures can be expressed in terms of the global temperature. 
However, the relationships between these partial temperatures are functions of the mole fractions (composition) and 
lead to new contributions (not considered in previous works) to the transport coefficients. The consequences of this 
new effect on transport are quite significant. 

The dependence of the transport coefficients on the full parameter space has been explored in the case of the shear 
viscosity 77. In addition, these theoretical predictions have been compared with the results obtained from a numerical 
solution of the Boltzmann equation by means of the DSMC method. The theory and simulation clearly show that 
in general, the influence of dissipation on 77 is quite important since there is a relevant dependence of the shear 
viscosity on the restitution coefficient. With respect to the accuracy of the Chapman-Enskog results we see that the 
first Sonine approximation exhibits an excellent agreement with the simulation data. This supports the idea that 
the Sonine polynomial approximation for granular fluids has an accuracy comparable to that for elastic collisions. 
Exceptions to this agreement are extreme mass ratios and strong dissipation. These discrepancies are due basically 
to the approximations introduced in applying the Chapman-Enskog, and more specifically in using the first Sonine 
approximation. 

Although the utility of a hydrodynamic description for granular media under rapid flow conditions has been recog- 
nized for many years, there are some doubts about the possibility of going from a kinetic theory to a hydrodynamic 
level of description by using a Chapman-Enskog expansion around the homogeneous cooling state. Given that the 
search for exact solutions of the Boltzmann equation is far beyond the present perspectives, an alternative to get some 
insight into the above question is to numerically solve the kinetic equation by using for instance the DSMC method 
and compare these results with the corresponding solution obtained by assuming the validity of a hydrodynamic 
description. The good agreement obtained here for the temperature ratio and the shear viscosity coefhcient shows 
clearly the direct verification of hydrodynamics and the quantitative predictions for transport coefficients from kinetic 
theory for states with small spatial gradients but including strong dissipation. 

One of the main limitations of the results obtained here from the Boltzmann equation is its restriction to the 
low-density regime. In this regime, the collisional transfer contributions to the fluxes are negligible and only their 
kinetic contributions are taken into account. Possible extension in both aspects, theory and simulation, of the present 
hydrodynamic description to higher densities can be carried out in the context of the revised Enskog theory. In this 
case, many of the phenomena appearing in dense granular fluids (such as spontaneous formation of dense clusters 
surrounded by regions of low-density j23|) could be studied. 

The instability of the HCS found in the previous section may also seem to represent a limitation on the validity 
of hydrodynamics. However, it should be noted that the Chapman-Enskog expansion is not an expansion about 
the HCS but rather its local form, parametrized by the exact hydrodynamic fields. These fields obey the nonlinear 
hydrodynamic equations which are not unstable, but rather they describe the correct subsequent evolution of the 
unstable linear equations. The linearized hydrodynamics is indeed limited to time scales short compared to the 
growth of the linear perturbations. Nevertheless, these equations describe correctly the nature of the onset of the 
instability. Since it is a long wavelength phenomenon, the instability can be suppressed completely for sufficiently 
small system sizes. 
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